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Addressing the spin question in gravitational- wave searches: 
Waveform templates for inspiralling compact binaries with nonprecessing spins 
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This paper presents a post-Newtonian (PN) template family of gravitational waveforms from inspiralling 
compact binaries with non-precessing spins, where the spin effects are described by a single "reduced-spin" 
parameter. This template family, which reparametrizes all the spin-dependent PN terms in terms of the leading- 
order (1.5PN) spin-orbit coupling term in an approximate way, has very high overlaps (fitting factor > 0.99) 
with non-precessing binaries with arbitrary mass ratios and spins. We also show that this template family 
is "effectual" for the detection of a significant fraction of generic spinning binaries in the comparable-mass 
regime (niT/nii < 10), providing an attractive and feasible way of searching for gravitational waves (GWs) 
from spinning low-mass binaries. We also show that the secular (non-oscillatory) spin-dependent effects in the 
phase evolution (which are taken into account by the non-precessing templates) are more important than the 
oscillatory effects of precession in the comparable-mass {nii ^ mi) regime. Hence the effectualness of non- 
spinning templates is particularly poor in this case, as compared to non-precessing-spin templates. For the case 
of binary neutron stars observable by Advanced LIGO, even moderate spins (L^. • S/rrr ^ 0.015 - 0.1) will 
cause considerable mismatches (~ 3% - 25%) with non-spinning templates. This is contrary to the expectation 
that neutron-star spins may not be relevant for GW detection. 



I. INTRODUCTION 

Coalescing compact binaries consisting of stellar- 
mass/intermediate-mass black holes and/or neutron stars 
are among the most promising sources of gravitational waves 
(GWs) for the interferometric GW detectors like LIGO, 
Virgo and GEO 600. Compact binary systems consisting of 
neutron stars or black holes can be produced in a variety 
of astrophysical scenarios, which can be broadly classified 
into two main classes: 1) Isolated binary evolution in which 
two massive stars constituting a binary undergo successive 
supernova explosions without disrupting the binary orbit [1] 
2) Dynamical formation scenarios in which two compact 
objects form a bound orbit due to dynamical interaction in 
dense stellar environments [2]. Once formed, the binary loses 
orbital energy and angular momentum through GW emission 
and starts to inspiral and finally the binary components merge 
with each other 

While spin measurements of stellar-mass black holes have 
indicated that many black holes may have very high spins 
(||S||/m^ ~ 0.2 - 0.98) [3], most of the observed neutron stars 
are found to be weakly spinning [4]. Additionally, for neutron 
stars, there is a theoretical upper limit (||S||/m^ ~ 0.7) on the 
spin rate beyond which the neutron star is gravitationally un- 
stable [5, 6]. In the case of compact binaries formed in isolated 
evolution, it is likely that the spins will be nearly aligned to the 
orbital angular momentum [7]. But the distiibution of the spin 
tilt angles (angle between the spin and orbital angular momen- 
tum) depends on the distribution of supernova kicks [8]. On 
the other hand, in the case of binaries formed in dynamical 
interactions, there is no prior reason to expect any spin align- 
ments. If the spins are misaligned with the orbital angular mo- 
mentum, the general relativistic spin-orbit and spin-spin cou- 
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pling will cause the spins to precess around the (nearly fixed) 
direction of the total angular momentum [9]. The complexity 
of the binary's dynamics will be encoded in the GW signals 
observed by a detector. 

The GW signals from compact binaries buiied in the noisy 
data of interferometric detectors are best extracted by employ- 
ing the technique of matched filtering, which involves cross 
coiTelating the data with theoretically calculated templates of 
the expected GW signals. Theoretical GW templates can be 
constructed by solving Einstein's equations using analytical 
approximation methods and/or numerical techniques. When 
the compact objects are well separated and are slowly mov- 
ing {vjc <K 1) gravitational waveforms can be computed using 
the post-Newtonian (PN) approximation to General Relativ- 
ity [10]. The PN approximation breaks down as the compact 
objects reach the ultra-relativistic regime, and an accurate de- 
scription of the merger process requires exact solutions of Ein- 
stein's equations which can only be obtained by large-scale 
numerical simulations [11-13]. 

Analytical GW templates (parameterized by the masses and 
spins) can be constructed by combining PN calculations with 
numerical-relativity simulations (see, e.g., [14-19]). In the 
case of "low-mass" {m\ + m2 < 12Mq [20-22]) binaries 
where the signal-to-noise ratio (SNR) is almost entirely con- 
tributed by the inspiral portion of the signal, it is sufficient 
to model only the inspiral accurately, where the PN approxi- 
mation holds. Although the expected form of the signals can 
be computed as a function of the source parameters, the pa- 
rameters of a signal that is buried in the noise is not known a 
priori. Thus the data is cross correlated with a "bank" of the- 
oretical templates corresponding to different (astrophysically 
plausible) values of physical parameters. A geometrical for- 
malism has been developed for "laying down" templates in the 
parameter space of compact binaries [23, 24]. Parameters of 
the templates are chosen in such a way that the loss of SNR due 
the mismatch between two neighboring templates is less than 
an acceptable value, while at the same time keeping the total 



number of templates employed in the search computationally 
tractable. 

However, templates for spinning binaries are characterized 
by a large number of parameters (two for the masses, six 
for the spins, two for the inclination and polarization angles), 
complicating placement of templates in the bank and consid- 
erably increasing the computational cost of searches. Further- 
more, many different spin configurations are known to be de- 
generate, making it unnecessary to employ templates corre- 
sponding to all parameters. 

The idea of constructing an effective template family pa- 
rameterized by a smaller number of parameters that has high 
enough overlaps with the expected signals was first proposed 
by Apostolatos [25], who introduced a modulational sinu- 
soidal term in the frequency-domain phase of the templates 
to capture the oscillatory effects of precession. However, it 
was found that this template family fails to capture the tar- 
get signals with sufficient efficiency [26]. Later, Buonanno, 
Chen and Vallisneri (BCV) [27] proposed a phenomenolog- 
ical template family that has better overlaps with the tar- 
get signals, by introducing modulation effects in both the 
frequency-domain- amplitude and phase of the templates. Sev- 
eral of these phenomenologically introduced parameters could 
be searched over analytically, thus reducing computational 
costs significantly. A search for GWs from spinning binaries 
was performed using the data from LIGO's third science run, 
employing BCV templates [28, 29]. Interestingly, the sensi- 
tivity of this search towards spinning binaries was found to 
be not any better than a search employing non-spinning tem- 
plates. Indeed, BCV have cautioned that the increased degrees 
of freedom associated with the detection statistic (due to the 
analytical maximization of several phenomenological param- 
eters) will also increase the false alarm rate. Later Van Den 
Broeck et al. [30] demonstrated that, while BCV templates 
have high overlaps with the target signals, the increased de- 
grees of freedom associated with the detection statistic and 
the lack of good methods for vetoing non-Gaussian detector 
glitches that mimic the expected signals ("signal-based ve- 
toes") resulted in producing a much higher false alarm rate 
in the presence of non-Gaussian noise, and hence negated 
the advantages. They argued that the standard non-spinning 
frequency-domain 3.5PN templates have higher detection ef- 
ficiency (at a given false alarm rate) towards spinning bina- 
ries than the current implementation of the phenomenological 
BCV template bank. 

BCV [27] also proposed a physical template family, which 
presumably will not suffer from this limitation. This template 
family, which assumes that only one compact object has sig- 
nificant spin, is parametrized by four parameters that need to 
be searched over (two masses, magnitude of the single spin, 
spin tilt angle). This template family was extensively studied 
by Pan et al. [3 1 ] and Buonanno et al. [32], where they demon- 
strated that the template family is effectual [33] in detecting 
single-spin as well as double-spin binaries using Initial detec- 
tors. They also demonstrated the possibility of reducing the 
number of spin parameters to one in certain regions in the pa- 
rameter space. This template family has been implemented in 
the LIGO- Virgo search pipelines [34-36]. Another enhance- 



ment of the non-spinning frequency-domain template family 
which extends the symmetric mass ratio (;; = mim2/in^) to 
unphysical values (77 > 0.25) has also been proposed as a de- 
tection template family for spinning binaries [37-39]. 

In this paper we propose a frequency-domain FN template 
family characterized by a single "reduced-spin" parameter 
(and two masses) for the case of binaries with non-precessing 
spins (spins aligned/anti-aligned with the orbital angular mo- 
mentum). This template family, which reparametrizes all 
the spin-dependent PN terms in terms of the leading-order 
(1.5PN) spin-orbit coupling term in an approximate way, has 
very high overlaps with non-precessing binaries with arbi- 
trary mass ratios and spins. We also show that this "reduced- 
spin" template family is able to capture a significant fraction 
of generic precessing binaries in the comparable-mass regime 
(q = m2lm\ < 10). This might provide an efficient and fea- 
sible way of searching for spinning low-mass binaries in the 
comparable-mass regime. This result also shows that the secu- 
lar (non-oscillatory) spin-dependent effects in the phase evolu- 
tion (which are taken into account by the non-precessing tem- 
plates) are more important than the oscillatory effects in the 
case of (nearly) equal-mass binaries even with generic spins. 

This paper is organized as follows: Section 11 summarizes 
the main findings of this paper as well as lists the limitations of 
this work. Section 111 provides a description of the PN wave- 
forms from precessing binaries computed in the adiabatic ap- 
proximation, which are assumed to be the "target signals" that 
we want to detect. Section IV is a brief introduction to the 
GW data analysis of inspiralling compact binaries. In Sec- 
tion V, we construct the frequency-domain PN template fam- 
ily parametrized by a reduced spin parameter and demonstrate 
the effectualness [33] of the template family in detecting bina- 
ries with non-precessing spins. In Section VI we investigate 
the effectualness of the template family in detecting precess- 
ing binaries. We use geometrical units throughout the rest of 
the paper: G - c - 1 . 



II. SUMMARY OF RESULTS, RECOMMENDATIONS, AND 
LIMITATIONS OF THIS WORK 

A brief summary of the main findings of this paper is given 
below: 

1 . For the case of binary neutron stars observable by Ad- 
vanced LIGO, even moderate spins (L^-S/m^ ^ 0.015- 
0.1) will cause significant mismatches (3% - 25%) with 
non-spinning templates. This is contrary to the expec- 
tation that spin effects may not be relevant for the de- 
tection of GWs from binary neutron stars. Secular spin- 
dependent effects in the phase evolution are more impor- 
tant than oscillatory effects in the case of comparable- 
mass {m\ ^ mi) binaries. Hence the effectualness 
of non-spinning templates is particularly poor in the 
comparable-mass regime, as compared to templates de- 
scribing non-precessing spins, which take into account 
the secular effects (see Figs. 9 and 1 1). 

2. In the case of binaries with non-precessing spins, it is 



possible to describe the dominant spin effects in terms 
of a single "reduced-spin" parameter (the combination 
of spins and mass ratio that describes the leading-order 
spin-orbit-coupling term/? [40]). The reduced-spin tem- 
plate family has very high fitting factor [41] towards 
non-precessing binaries with arbitrary spins and mass 
ratios (see Fig. 5), while the non-spinning template fam- 
ily can cause significant loss of SNR in certain regions 
of parameter space, e.g., spins aligned to the angular 
momentum (see Fig. 4). Note that binaries with spins 
nearly aligned to the angular momentum are expected 
from isolated binary evolution [7]. 

3. Assuming a target population of binaries with uni- 
formly distributed spin magnitudes (||S||/ot^) in the in- 
terval (0, 0.98) for black holes and (0, 0.7) for neu- 
tron stars, and isotropically distributed spin angles and 
inclination-polarization angles, we show that almost the 
entire population of equal-mass binaries with generic 
spins can be detected using the reduced-spin templates 
with very high fitting factors. On the other hand, only 
~ 51% (52%) 57% of the equal-mass population with 
total mass 2.^Mq(12Mq)20Mq produces the same fit- 
ting factor with non-spinning PN templates. Even if we 
restrict the neutron star spins to (0, 0.3), this result does 
not change significantly (see Fig. 9). 

4. The effectualness of the reduced-spin templates de- 
creases with increasing mass ratio (m2/»Ji » 1), 
since for binaries with significantly unequal masses, os- 
cillatory effects of precession become crucial. Still, 
the reduced-spin template family performs considerably 
better than non-spinning templates: ~ 60% (45%) 37% 
of the population with milnii - 3.5 (7.1) 13.3 produces 
fitting factor > 0.97 with reduced spin templates, while 
~ 44% (36%) 32% produces the same fitting factor with 
non-spinning templates (see Fig. 9). 

5. In the case of unequal-mass binaries, most of the bina- 
ries producing fitting factor < 0.9 with the reduced spin 
templates are significantly tilted with respect to the de- 
tector and have spins of the more massive object highly 
non-aligned with the orbital angular momentum. In the 
case of GW searches where astrophysical priors restrict 
the inclination angles or spin tilt angles of the target 
population to small values (e.g., binary mergers produc- 
ing short-hard gamma ray bursts observable to us), it 
might suffice to consider only the secular spin effects 
(see Fig 10). 

6. For highly unequal-mass binaries, spin effects are al- 
most entirely determined by the spin of the more mas- 
sive compact object [89], and hence a template family 
describing precessional effects assuming only one spin- 
ning compact object, such as the physical template fam- 
ily proposed by BCV [27, 31, 34, 35] should be able to 
model these binaries accurately. This suggests one natu- 
ral way of splitting the parameter space in GW searches: 
Binaries in the comparable-mass regime (where pre- 
cession effects are not significant) can be detected em- 



ploying templates described by a single spin parameter, 
while the detection of binaries in the high mass ratio 
regime (where precession effects are significant) will re- 
quire templates described by two spin parameters. 

7. The observation that binaries in the comparable -mass 
regime can be detected employing a template family de- 
scribed by a single spin parameter also suggests that it 
will be hard to estimate the other spin parameters of the 
binary accurately employing restricted PN templates. 
This underlines the need of including all the known 
physical information (from higher harmonics, merger- 
ringdown etc.) in the parameter-estimation pipelines 
which might help to disentangle the correlations be- 
tween different spin components. 

A note on the limitations of this work: We consider only 
the dominant harmonic of the gravitational waveforms from 
the inspiral stage of the binary coalescence, and the amplitude 
of the waveforms is computed only in the leading order ("re- 
stricted" PN approximation). In the detection problem, higher 
harmonics are believed to be not important for the case of low- 
mass binaries in the comparable mass regime. But for more 
massive and highly unequal-mass binaries, the contribution 
from higher harmonics can be significant - especially if the 
binary is tilted with respect to the detector [42, 43]. Also, the 
effect of post-inspiral stages is negligible in the case of low- 
mass binaries (m < 12Mq), but can not be neglected in the 
case of high-mass binaries. Still, we consider binaries with m 
as large as 2OM0 in this paper. The rationale for this choice is 
that in this paper we are investigating only the spin-effects in 
the GW detection problem. Since we are neglecting the effects 
of post-inspiral stages, we evaluate the performance of differ- 
ent template families only in the frequency range / < /isco, 
where /isco is the frequency of the innermost stable circular 
orbit (ISCO) in the Schwarzschild geometry. Note that, for 
the case of spinning binaries, the "actual" ISCO frequency de- 
pends on the spin; we ignore this. 

Additionally, in this paper, we quantify the effectualness of 
the template family only in terms of the fitting factor, which 
quantifies the loss of SNR due to the mismatch of the signal 
and template in Gaussian noise. Since actual detector data is 
not perfectly Gaussian, additional signal-based vetoes are used 
in the detection statistic, such as the "chi-square veto" [44], 
whose effects are not considered here. But, note that, un- 
like the case of the BCV phenomenological template family, 
signal-based vetoes can be readily implemented for the case 
of the reduced-spin template family. But we note that a com- 
plete characterization of the effectualness of a template family 
requires detailed studies using actual detector data. This is be- 
yond the scope of this paper. 



ni. POST-NEWTONIAN WAVEFORMS IN THE 
ADIABATIC APPROXIMATION 

In this discussion, we will only consider binaries in quasi- 
circular [90] orbits. Indeed, inspiralling compact binaries 



are expected to lose their eccentricity before their GW sig- 
nals reach the sensitive frequency band of ground-based de- 
tectors [45, 46]. In the early stages of the evolution of the 
binary due to GW emission, the change in orbital frequency 
is much smaller than the orbital frequency itself. During this 
adiabatic inspiral, the loss of the specific orbital binding en- 
ergy E{v) (binding energy per unit mass) is related to the en- 
ergy flux of gravitational radiation Tiv) in the following way: 
mdE{u)/dt = -Tiv). This "energy balance" argument pro- 
vides the following coupled ordinary differential equations 
from which orbital phase evolution ip{t) of the binary can be 
computed 



dip 
~dt 



V 

m 



du _ Tip) 
dt mE'(v) 



(3.1) 



Above, y is a velocity parameter which is related to the orbital 
frequency w by y = (mw)'^-' where m = m\ + m^ is the total 



mass of the binary, and E'(v) = dE(v)/dv. The energy and 
flux functions can be computed as PN expansions in terms of 
the small parameter v. Currently the energy function has been 
computed up to 3PN order (v^) [42, 47-60] and flux function 
up to 3.5PN order (v^) [42, 53, 54, 56, 61-63]; but the spin 
effects have been computed only up to 2PN and 2.5PN orders, 
respectively [91]. 



The energy balance equation can be solved in a number of 
different ways (see [22] for a recent overview) which are per- 
turbatively equivalent at the corresponding PN order. But, 
since the PN expansion of the energy and flux functions is 
known up to only a limited PN order, the actual results of these 
different methods can be somewhat different. In this paper, we 
construct a new approximant to the GW phasing computed by 
re-expanding the rational function E'{v)l'F{v) and truncating it 
at 3.5PN order (spin terms are considered only up to 2.5PN); 
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Above, 7/ = m\m2lm is the symmetric mass ratio, 6 = 
{nil - m2)/m the asymmetric mass ratio, je the Euler's con- 
stant, and La, the unit vector along the Newtonian orbital an- 
gular momentum La, = m^rj/u h^. The spin variables Xs and 
Xa are related to the (dimensionless) spin vectors of the binary 



as: 



Xs = (Xi + X2)/2, Xa = (Xi - X2)/2, 



(3.3) 



where Xi - Si/m^, S, being the spin angular momentum of the 
/th compact object. 



Equation (3.2) can be plugged back to Eq. (3.1) to get an 
explicit expression of the orbital phase in terms of v. 
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where ipQ is a certain reference phase. 

Thus the phasing formula can be reduced to one differential 
equation describing the evolution of the orbital frequency, and 
one explicit expression of the orbital phase: 
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We call this particular way of solving the phasing formula the 
"TaylorT5" approximant [92]. The choice of this particular 
approximant as the target signal family is motivated by the 
following reason: Since the spin-dependent terms in the PN 
expansion of the energy and flux functions are available only 
up to a rather low 2.5PN order, the different approximants give 
somewhat different results. We want to isolate this issue from 
the issue of the effect of spin-precession in the target signals. 
Thus, we construct the target waveforms in such a way that 
they are as close to the non-precessing frequency-domain tem- 
plate family as possible in the limit of non-precessing spins. 
Since the frequency-domain template family used in this pa- 
per ("TaylorF2" approximant [22]) is constructed based on a 
re-expansion of E'{u)/f{u), we choose to construct the time- 
domain target waveforms also based on this re-expansion. 
Note that the effectualness of the template family (although 
weakly) depends on the particular approximation used in the 
construction of the target and template waveforms. This is an 
indication of the level of truncation error in the PN expansion, 
and points to the need of computing the higher PN order spin 
terms. This is being explored in an ongoing work [64]. Also 
note that the "TaylorT5" approximant has an additional advan- 
tage (over TaylorTl and TaylorT4) that only one differential 
equation needs to be solved numerically; the orbital phase is 
computed as an explicit expansion in u. 

If the spin vectors are misaligned with the orbital angular 
momentum, the spin-orbit and spin-spin coupling cause the 
spins and orbital angular momentum to precess around the 
nearly constant direction of the total angular momentum J, 
constantly changing the angle between the spins and angular 
momentum [9]. The evolution equations for the orbital angu- 
lar momentum and spins, including the next-to-leading-order 
spin-orbit terms, are given by [27, 53]: 
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Above, q = m2lm\ is the mass ratio. The instantaneous pre- 
cession frequency of the individual spins is ||il,||. 

The orbital frequency, spins and orbital angular momentum 
can be evolved by solving the differential equations Eqs. (3.5), 
(3.6) and (3.7). In order to perform the evolution, we adopt 
the coordinate system proposed by Finn and Chernoff [65], 
and subsequently used by several authors. The z-axis of this 
(fixed) source coordinate system is determined by the initial 
total angular momentum vector J ^ L^r -H Si -H S2 and the x- 
axis is chosen in such a way that the detector lies in the x- z 
plane. 

In the absence of precession, the phase evolution of the 
(dominant harmonic) GWs is twice the orbital phase if(t). But 
the precession of the orbital plane introduces additional mod- 
ulation in the orbital phase. If we define Q>{t) as the orbital 
phase with respect to the line of ascending nodes (the point at 
which the orbit crosses the x - y plane from below) then the 
the phase evolution of the (dominant harmonic) GWs is given 
by 20(f), where 0(f) is given by 



JO 

dt 



da 
dt 



■ cos 7, 



(3.9) 



where a and / are the angles describing the evolution of the 
orbital angular momentum vector L^ in the Finn-Chemoff co- 
ordinate system: 



a = arctan(LAfy/LAfv), i = arccos(LAij.). 



(3.10) 



Computation of gravitational waveforms in the detector frame 
follows the description of BCV [27] (Section IIC). In the re- 
stricted PN approximation, the resulting gravitational wave- 
form observed at the detector can be written as: 

h(t) = CQ{t) cos2 [0(f) + Oo] + 5Q(f) sin2 [0(f) + O,,] , 

(3.11) 



where 

Cait) = -^t;Mc+(0^+ + Cx(0^x] , 

Sait) = -^v" [S4t)F^+SAt)F^] , (3.12) 

and (to is a constant that depends on the initial configuration of 
the binary. Above, D is the luminosity distance to the binary, 
while F+ and Fx are the antenna pattern functions: 



F^ 



( 1 + cos 0) COS 20 cos 2)/' - cos 6 sin 20 sin 2i/r , 



Fx - -(1 + cos 0) cos 20 sin 2i/' + cos sin 20 cos 20, 

(3.13) 

where 9 and are the polar and azimuth angles specifying 
the position of the binary in the detector frame, and is the 
polarization angle. Also, in Eq. (3.12), 
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where is the angle between the line of sight from the detector 
to the source and the z-axis of the coordinate system (direction 
of initial J). Note that, while is a constant, the angles a and 
/ vary throughout the evolution of the binary (see Fig. 1). 

The expression for the observed GW signal in a detector 
given in Eq. (3.11) can be rewritten as 



h(t) = A(t) cos imt) + 'Pol . 



(3.15) 



where 



A(t) = ^Cl(t)+Sl(t), 



^'{t) = 20(0 - arctan 



SqU) 



CqH) 



(3.16) 



and 'i'o is a constant that depends on the initial configuration of 
the binary. Note that S git) and Cgit) are oscillatory in nature. 
These expressions show that the precession of the orbital plane 
is causing amplitude- and phase modulation in a GW signal 
h(t) observed by a fixed detector 



A. Initial conditions and truncation conditions 

In order to simulate a GW signal observed in a detector, we 
have to specify, apart from the masses and spins, the initial 




► ey 



FIG. 1: Source frame in the Finn-Chemoff convention. The z-axis 
points to the direction of the initial total angular momentum vector J, 
the detector lies in the x- z plane, and the angles a and i describe the 
evolution of the orbital angular momentum vector L^v in the source 
frame. 



orbital phase ipo, two angles and describing the initial ori- 
entation of the binary in the detector frame, and two angles 
and describing the sky-location of the binary in the detector 
frame. In this paper, we will assume that the binary is located 
optimally with respect to the detector: = = 0. Also, we 
set the initial phase tpo = 0. The evolution of the binary is 
started at the specified low-frequency (20 Hz) and truncated 
when E'iv)/r(v) = [93]. 

Spins are specified in a coordinate system whose z-axis is 
determined by the initial Newtonian orbital angular momen- 
tum Lfti. After computing J ^ L^ + Si -H S2, we rotate 81,82 
and Lat such that initial J is pointed along the z-axis. 

It is often convenient to describe the spin vectors x = 
iX-fXy^Xz) in terms of a spin magnitude ||x|| and two spin an- 
gles 6s and 0^ : 

Xx = 11x11 sine's cos 05, 
Xy = 11x11 sine's sin 0s, 
X,. = 11x11 cos 0s, (3.17) 

where 0s is referred to as the spin tilt angle. 



IV. DETECTION OF BINARY INSPIRAL SIGNALS 

The GW signal h{t; A) observed in a detector is a function 
of the set of physical parameters A of the binary, such as the 
component masses and the spins. Detecting the signal requires 
analyzing (noisy) interferometric data. We assume that the 
detector noise n{t) follows a zero-mean Gaussian distribution, 
characterized by its (one-sided) power spectral density (PSD) 
S hif)- We also assume the noise to be additive. This implies 



that when a signal is present in the data d(t), then 

d(t) = h(t; X) + n(t) . (4.1) 

Under the above assumptions about the characteristics of de- 
tector noise, the Neyman-Pearson criterion [66] leads to an 
optimal search statistic, which when maximized over the over- 
all amplitude of the signal, is the cross-correlation of the data 
with a normalized template. 



p = {h(X), d) 



(4.2) 



where the normalized template is h{f) = h{f)l y/ih, h). The 
angular bracket denotes the following noise-weighted inner 
product, called the overlap: 



{a, b) = 451 



-'/to, 






(4.3) 



where 5(/) and b(f) are the Fourier transforms of a{f) and b{f), 
respectively. Also, /low is a low-frequency cutoff of the detec- 
tor noise and /cut is an upper cutoff frequency where the tem- 
plates cease to describe the true signal with sufficient accuracy 
(e.g., due to the PN approximation breaking down). The opti- 
mal SNR of the filter is given by 



Pop, = </2(A), /2(A)) 



1/2 



(4.4) 



In a "blind" search in detector data, where none of the bi- 
nary's parameters are known a priori, the search for a GW 
signal requires maximizing p over a "bank" of templates cor- 
responding to different values of A. Apart from the physical 
parameters A, the waveform also depends on the (unknown) 
time of arrival to of the signal at the detector, and the corre- 
sponding phase *Po. Maximization over *Po is effected by using 
two orthogonal templates for each combination of the physi- 
cal parameters [67], and the maximization over fo is attained 
efficiently with the help of the Fast Fourier Transform (FFT) 
algorithms. 

The inner product in Eq.(4.3) between the template x(A) 
(described by the set of parameters A) and the target signal 
h{\) maximized over 'i'o and fo is called the match: 



match - max4/|,,,„ {xi^a, fo. A), h{\)). 



(4.5) 



The fraction of optimal SNR retrieved by a sub-optimal tem- 
plate bank is given by the match maximized over all the tem- 
plate parameters, called ih& fitting factor [41]: 



FF = maxvp„,,„,A {xQVq, fo. A), h{\)). 



(4.6) 



Template families with high fitting factors (FF > 0.97) are con- 
sidered to be "effectual" for signal detection, while templates 
with high values of match (match > 0.97) are considered to 
be both effectual in detection and "faithful" in estimating the 
parameters [33]. 



A. Computing overlaps 

Computing the inner product in Eq.(4.3) requires modeling 
the PSD of the detector noise. In this paper, we use a fit to 



the expected PSD of the "zero-detuning, high power" (zero- 
detuning of the signal recycling mirror, with full laser power) 
configuration [68] of Advanced LIGO [69]: 

S,,{f) = 10"'*'*(0.01 52 ;c"'*H- 0.2935 /^"^-H 2.795 1jc^^2 

- 6.5080 x^''' + 17.7622) , (4.7) 

where x - //245.4. The low-frequency cutoff (due to the 
elevated seismic noise) is assumed to be 20 Hz. 

The upper cutoff frequency of the integral in Eq.(4.3) is 
chosen as the frequency of the innermost stable circular orbit 
(ISCO) of a test particle around a Schwarzschild black hole 

/isco = ^isco/'^-''^'"-'' where uisco = 1/ V6. Note that, for the 
case of comparable-mass binaries, ISCO is a poorly defined 
quantity (see, e.g. [70]), and the validity of the PN approxi- 
mation can be only tested by comparing against fully general- 
relativistic numerical simulations (see [71-77] for some work 
in this direction). Nevertheless, /isco can be treated as a con- 
venient cutoff frequency, especially since we will be comput- 
ing overlaps with two different template families (spinning and 
non-spinning), and a fixed upper cutoff frequency will pro- 
vide a fair comparison between the two. Indeed, in the case of 
spins (nearly) anti-aligned with the orbital angular momentum, 
the "actual" ISCO frequency could be significantly lower than 
/isco, and the time-domain waveforms discussed in Sec. Ill 
will terminate before /sco- This will reduce the overlaps of 
the templates with the target signals. But, since different tem- 
plate families will be affected in the same way, we neglect this 
effect. But we note that, in a actual search it might be more 
appropriate to use a spin-dependent cutoff frequency. Fourier 
transform of the time-domain target signals is computed with 
the help of the FFTW library [78]. In order to avoid artifacts as- 
sociated with the abrupt start and stop of the waveform in the 
time domain, a tapering window [79] is applied prior to com- 
puting the FFT. Maximization of the overlaps over the physi- 
cal parameters (masses and spins) is performed with the aid of 
the Nelder-Mead downhill simplex algorithm amoeba [80, 81] 
with the implementation described in the Section III B of [82], 
except that the mass parameters used in the maximization pro- 
cedure are the chirp mass M = m 77^^^ and 77. 



V. FOURIER DOMAIN TEMPLATES FOR INSPIRALLING 
BINARIES WITH NON-PRECESSING SPINS 

In the case of binaries with spins aligned/anti-aligned with 
the orbital angular momentum, it can be seen from Eq. (3.7) 
that dSj/dt = 0, and hence the spins and the angular momen- 
tum do not precess. This will considerably simplify the com- 
plexity of the GW signal, and all the spin effects can be de- 
scribed by two parameters Xs = Xs ■ i^N = Of 1 + ^2)/2 and 
Xa = Xa' i^N = (Xi ~X2)/'2; which remain constant throughout 
the evolution. 

Since the cross correlation between the data and the tem- 
plate is most efficiently computed in the Fourier domain by 
using the FFT, waveform templates in the Fourier domain are 
computationally cheaper In the case of binaries with non- 
precessing spins, Fourier transform of the time-domain tem- 



plates can be computed analytically using the stationary phase 
approximation [83], which can be written as 



Above, 



h{f;XuX2)=A(f)e 



Mn-c'J^} 



-i[2;r//o+4'o+'i'(/)-;r/4] 



dF(vf) 



dt 



-1/2 



(5.1) 



(5.2) 



v\l{nm) is the instantaneous 



where Vf - (mn/y^^, F(uf) 
GW frequency (of the dominant harmonic) evaluated at the 
stationary point Vf, C is a numerical constant that depends on 
the relative position and inclination of the binary with respect 
to the detector (C = 1 for optimally located and oriented bina- 
ries), to is the time of arrival of the signal at the detector, 'i'o 
the corresponding phase, and [40, 42] 
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back to Eq.(5.2) to get a closed form expression of A(/). 



denote the leading-order spin-orbit coupling, leading-order 
spin-spin coupling, and next-to-leading-order spin-orbit cou- 
pling, respectively. 

The time-derivative of the GW frequency appearing in 
Eq.(5.2) can be written as: 



(5.5) 



where we have used the energy balance equation indE/dt - 
—Tiv), and the definition E'iv) = dE/dv. This can be plugged 
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Using the re-expansion of E'(v)l'F{v) given in Eq.(3.2) we get 
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/ 502429 907// 



\ 16128 192 ) ^""^\ 48 



5rf 73921// 502429\ 



Xs- 



2016 "^ 16128 / 

(5.8) 

Note that the higher order terms in A(/) are not due to the 
PN corrections to the (time-domain) amplitude of the wave- 
form, but come from the higher order corrections to dF/dt. 
The standard practice followed in the literature is to truncate 
the expansion of dF/dt in the leading order But we found that 
including the higher order corrections to dF/dt gives a better 
agreement with the exact (numerical) Fourier transform of the 
time-domain signals. 

It is well known that the leading order spin term y6 (spin- 
orbit coupling) appearing at 1 .5PN order in the amplitude and 
phase can be represented by a single parameter (see e.g. [40]), 
which we call the "reduced spin" parameter: 



X =Xs + Sxa - YY^ATs 



(5.9) 



4) 
' Note that this is different from the "effective-spin" Sefi = 

(Xs + ^Xa - TlXsl"^) w^, introduced by Damour [84] to describe 
the leading order spin-orbit term in the conservative dynam- 
ics of the binary. In some sense, the reduced-spin parameter ;t' 
describes both the conservative and radiative dynamics of the 
binary as encoded in the GW signal. 

The spin-dependent terms y6, cr, y and e defined in Eqs.(5. 10) 
and (5.8) can be written in terms of the reduced-spin parameter 
^as: 

P = 113;r/12, 
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FIG. 2: Relative significance of the spin terms neglected in the 
"reduced-spin" templates, as a function of the symmetric mass ra- 
tio r]. The plot shows the ratios of the coefficients of the spin terms 
Xs andxci (which are neglected in "reduced- spin" templates) with the 
coefficients of ;(' (which are included). It can be seen that the rel- 
ative contribution of the neglected terms is small (< 10%), and the 
dominant spin effects are described by the single parameter x- See 
Eq.(5. 10) for the definition of cr,, 7, and e,. 
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Figure 2 illustrates the fact that the coefficients (0-1,0-2, 71, ei) 
of Xa and Xs ^e significantly smaller than the coefficients 
(0-0,70,^0) of;^f. This gives the indication that the dominant 
spin effects in the waveform could be captured by just using 
the reduced spin parameter x- A template family solely de- 
scribed by the reduced spin parameter x, apart from m and rj, 
can be constructed by setting 0-1 = 0-2 = 71 = ei = in 
Eqs.(5.1)-(5.11). i.e., 

Mf;x) = Kf;XuX2), with o-i = 0-2 = 71 = fi = O (5.12) 
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FIG. 3: Contribution to the phase 4'(/) at ISCO from the spin terms 
neglected in the "reduced-spin" template family. The plots show the 
contours of lA^Piscol (see Eq. (5.13) for definition) as a function of 
the spin parameters xi ^.ndxi of the two compact objects (xi being 
the spin of the more massive object), for the case of six different mass 
ratios. Darker shades correspond to lower values. 



It is possible to get a rough idea of the effect of the neglected 
spin terms 0-1,0-2,71, ei at different regions in the parameter 
space by looking at their contribution to the Fourier domain 
phase at the ISCO frequency: 



A^F 



-3 
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^0(o-lXa+(^2XsXa) 



fisco 



128// 
+ yiXa{'^ + 3 ln(!;isco)) 



(5.13) 



In Fig. 3, we plot the contours of lA^FiscoL which suggests 
that, for the equal-mass (77 - 0.25) case there is little loss in ne- 
glecting the additional spin terms 0-1,0-2,71 and e\ over most 
regions of the spin parameter space; but the differences get 
larger for the case of highly unequal masses. Nevertheless, it 
is likely that a template family described by one reduced spin 
parameter x will be able to capture the more general family 



10 



of signals (described by two spin parameters Xs ^ndxa) at the 
cost of a parameter bias in x- 

Figure 4 shows the match (see Eq. 4.5 for definition) of the 
"reduced-spin" template family h{f',x) with h(f',Xi^X2) for 
different values of the mass parameters m\ and m2 and spin 
parameters x\ ^i^'^ X2- It can be seen that the plots follow the 
same trends shown by Fig. 3 - the "reduced-spin" template 
family produces very high matches (> 0.99) either when the 
masses are equal or when the spins are equal. On the other 
hand, for the case of highly unequal masses and spins, the 
match can be as low as 0.8. Figure 5 plots the fitting fac- 
tor (see Eq. 4.6 for definition) of the reduced spin template 
family. Fitting factor is greater than 0.999 over the entire pa- 
rameter space, including extreme spins. This suggests that 
the reduced spin template family can be used in the detection 
of non-precessing-spin binaries with arbitrary mass ratios and 
spins causing no appreciable loss ofSNR. 

Figure 6 plots the fitting factor of the non-spinning PN 
template family, constructed by setting (3 - cr = y = 
in Eqs. (5.1)-(5.10), in detecting non-precessing signals 
h{f;X\^X2), which demonstrates the effect of neglecting the 
spin terms. The biggest eff'ect is in the case of binary neu- 
tron stars (oti - m2 - I.4M0). Since the secular spin effects 
accumulate over the large number of GW cycles present in 
the detector band, even moderate spins can cause appreciable 
mismatches in the case of binary-neutron-star signals in Ad- 
vanced LIGO. Note that these eff'ects are fully captured by the 
reduced spin templates. The non-spinning template family is 
least efficient towards binaries with the spin of the more mas- 
sive object aligned (as opposed to anti-aligned) to the orbital 
angular momentum. 

Figure 7 plots the mismatches (1 -fitting factor) of the non- 
spinning template family with non-precessing signals in the 
case of binaries with spins aligned to the angular momentum, 
as function of the spin of the target binary. Mismatches are 
computed for binaries with different component masses. For 
simplicity, the binary components are assumed to have equal 
spins (x\ - Xi)- It can be seen that, for the case of binary neu- 
tron stars {mi < 2Mq), even moderate spins (xj ^ 0.015 - 0.1) 
can cause appreciable mismatches (3% - 25%). This is con- 
trary to the prevailing assumption that the neutron star spins 
may not be relevant for GW detection. 

Finally, let us note that the "reduced-spin" template family 
proposed in this paper is different from the frequency-domain 
template family considering only the leading order spin-orbit 
term /3. In reduced-spin templates, we are re-expressing the 
higher order spin terms (appearing at 2PN and 2.5PN order in 
amplitude and phase) in terms of the leading-order spin term in 
an approximate way, without increasing the dimensionality of 
the template bank. Indeed, the main reason for the improved 
eff'ectualness of reduced-spin templates (as compared to non- 
spinning templates) is the inclusion of the leading order spin- 
orbit term. Nevertheless, inclusion of the higher order spin 
terms does provide some additional increase in effectualness. 
As an illustration, we plot in Figure 8, the fitting factor of the 
frequency-domain template family considering only the lead- 
ing order spin term (but considering non-spinning terms up to 
3.5PN). A comparison with Figure 5 will demonstrate the ad- 



vantage of the reduced spin template family over the template 
family considering only the leading order spin effects. 



VI. EFFECTUALNESS OF THE REDUCED-SPIN 

TEMPLATE FAMILY IN DETECTING GENERIC 

PRECESSING BINARIES 



TABLE I: Parameters used for the Monte-Carlo simulations of pre- 
cessing PN binaries. 



Parameter 


Value 


||X,-|| (m, > 2Me) 


uniform(0, 0.98) 


||X,-|| (m, < 2Mo) 


uniform(0, 0.7) 


0s 


uniform(0, n) 


<ps 


uniform(0, 2k) 


®,ifr 


uniform(0, n) 


0,(p 





A'sims 


~ 2500 — 5000 



In this section, we investigate the effectualness of the 
reduced-spin template family in detecting generic precessing 
PN binaries. There is suggestive indication that we will be able 
to model the most important spin effects in the orbital phase 
ip{u) using the reduced spin parameter ;k' = (1 - 7677/113) Xs ■ 
Lw + SXa • L^r — the leading-order spin orbit coupling in 
if{u) and dv/dt is completely described by this parameter (see 
Eqs. (3.2), (3.4) and (3.5)). Indeed, for the case of precessing 
binaries, the angles between the spin and orbital angular mo- 
mentum vectors change during the evolution, and hence the 
value of ;^'. But it can be shown that, over a significant fraction 
of the parameter space that we are interested in, ;t' is a slowly 
evolving quantity. 

Let us consider the precession equations (3.7) and (3.8). 
Considering only the leading order terms, the change in;^ can 
be computed as: 



dx 
dt 



39 6 
226 m^ T] 



u^ S2 . (S, X L^) 



(6.1) 



It is evident that x (in the leading order) remains a constant 
throughout the evolution either when one of the spins are zero 
(S; — > 0) or when the spins are parallel to each other or to the 
orbital angular momentum (S2 ■ (S 1 x L^^) —^ 0), or when the 
masses are equal (S — > 0). Additionally, the change in x is 
of the order of v^, and hence smaller than the change in the 
spin vectors, which is of the order of u^ (see Eqs. (3.7) and 
(3.8)) [94]. 

This observation suggests that the reduced-spin template 
family described in Section V might be able to detect a large 
fraction of precessing binaries in the comparable mass regime. 
We expect the effectualness of the template family to deterio- 
rate for the case of highly unequal masses. Also note that the 
precession of the orbital plane (described by the normal vec- 
tor Lat) for a given spin configuration is inversely proportional 
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FIG. 4: Match of the reduced-spin template family h{f\x) with non-precessing-spin signals h{f\ Xi'X2) in Advanced LIGO noise spectrum. 
Horizontal axes report the spin xt of the lighter object and vertical axes report the spin xi of the more massive object. The title of each plot 
reports the component masses (mi , ma) in units of Mq. 
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FIG. 5: Fitting factor (match maximized over M, 77 and x) of the reduced-spin template family h(f;x) with non-precessing-spin signals 
h(f',Xi'X2)- See Fig. 4 for a full description. 



to the symmetric mass ratio, as seen in Eq. (3.6) by taking 
the leading order terms, which also suggests that precessional 
effects (of the orbital plane) are minimal in the case of equal- 
mass binaries. Indeed, the GW signal observed by a detector 
contains additional modulations due to the precession of the 
orbital plane, as described by Eqs. (3.9) and (3.16), which de- 
pend on the orientation angles and ifr of the binary apart from 
the spin magnitudes, spin orientations and the masses. 

In order to quantify the effectualness of the reduced-spin 
template family in detecting precessing binaries, we perform a 
Monte-Carlo simulation where we generate generic spinning 
binaries with random spins and compute the fitting factor of 
the reduced spin template family with these target signals. For 
binaries with m > 2Mq, spin magnitudes ||x,|| are randomly 
drawn from a uniform distribution in the interval (0, 0.98), and 
for binaries with m < 2Mq spin magnitudes are distributed in 
the interval (0,0.7). The spin angles 0s and (f)s are uniformly 



distributed in the intervals (0, n) and (0, 2n), respectively. Ad- 
ditionally, the inclination of the initial total angular momen- 
tum vector J with respect to the line of sight from the detector, 
and the polarization angle i// are also uniformly distributed in 
the interval (0, tt). In the simulations, the binaries are assumed 
to be optimally located in the sky {0 - (p -0). Around 2500 — 
5000 simulations were performed for each mass configuration. 
The simulation parameters are summarized in Table I. 

Figure 9 presents cumulative histograms showing the frac- 
tion of binaries producing fitting factors > FF with the 
reduced-spin PN templates (solid lines) and non-spinning PN 
templates (dashed lines). These plots show that almost the 
entire population of equal-mass binaries can be detected us- 
ing the reduced spin templates with fitting factors > 0.97 [95]. 
On the other hand, only ~ 51% (52%) 57% of the equal-mass 
population with total mass 2.8M0 (12Mq) 20Mq produces the 
same fitting factor with non-spinning PN templates, demon- 
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FIG. 7: Mismatch (1-FF) of the non-spinning PN template family 
with non-precessing-spin signals with equal spins. The horizontal 
axis reports the spins of the target binary (xi = Xi)- Component 
masses of the target binaries (in units of M©) are shown in the leg- 
end. 3% and 10% mismatch values are indicated by horizontal dashed 
(black) lines. 



strating the importance of taking into account the secular spin- 
dependent phase evolution of the binary. The loss of non- 
spinning templates is the highest for the case of the binary 
neutron star system (IAMq, IAMq), contradicting the expec- 
tation that spin effects may not be significant for the search 
for binary neutron stars. The effectualness of the non-spinning 
templates gradually improves with increase in total mass, ow- 
ing to the smaller number of GW cycles present in the detector 
band. As expected, the effectualness of the reduced-spin tem- 
plates decreases with increasing mass ratio ^ (// <K 0.25), since 
for binaries with significantly unequal masses, oscillatory ef- 



fects of precession become important. Still, the reduced- 
spin template family performs considerably better than non- 
spinning templates: ~ 60% (45%) 37% of the population with 
m2//«i = 3.5 (7.1) 13.3 produces fitting factor > 0.97 with re- 
duced spin templates, while ~ 44% (36%) 32% produces the 
same fitting factor with non-spinning templates. For highly 
unequal-mass binaries, precessional effects are almost entirely 
determined by the spin of the larger compact object, and hence 
a template family describing precessional effects assuming 
only one spinning compact object, such as the physical tem- 
plate family proposed by BCV [27, 31, 34, 35] should be able 
to model these binaries accurately. 

It is worthwhile to identify regions of the parameter space 
where the reduced spin template family is not effectual so that 
a more complex template family could be used in those re- 
gions. Figure 10 shows scatter plots of the fitting factor of 
the reduced spin templates as a function of various parameters 
describing the initial configuration of the precessing binaries, 
while Fig. 1 1 shows the same for the case of non-spinning tem- 
plates. The main conclusions one can draw from these figures 
are: 

1. In the case of equal-mass binaries, all the binaries pro- 
ducing fitting factor < 0.9 with the reduced spin tem- 
plates correspond to the initial orientations of ^ 7r/2 
and i// ^ n/A, producing little or no observed signal in 
the detector. 

2. In the case of unequal-mass binaries, most of the bi- 
naries producing fitting factor < 0.9 with the reduced 
spin templates are significantly tilted with respect to the 
detector (0 » 0) and have spins of the more massive 
object highly non-aligned with the orbital angular mo- 
mentum (1 1:^2 1 1 sin 0s 2 s> 0). If there are strong astro- 
physical priors restricting the inclination angles or spin 
tilt angles of the target population to small values, the 
effectualness of this template family will be better than 
that presented in Fig. 9. A possible scenario where this 
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FIG. 9: Fraction of precessing PN binaries producing fitting factor 
> FF with the reduced spin PN templates (thick red solid lines) and 
non-spinning PN templates (thick red dashed lines). The total mass 
m/Mo and mass ratio q are shown in the titles. Thin black lines plot 
the same except that the spin distribution of neutron stars is restricted 
to the interval (0, 0.3). 



could apply is the GW followups of short-hard gamma 
ray bursts, where the opening angle of the jet is con- 
strained to a few tens of degrees (see, e.g., [85]) [96]. 



3. For the case of non-spinning templates (Fig. 11), the 
initial value of the reduced spin parameter x = {^ - 
76?;)/ 113 X.5 ■ Liv -H ^Xa ■ Lai clearly separates the popula- 
tion of binaries producing poor fitting factor from those 
producing high fitting factor. In particular, in the case of 
equal-mass binaries, those with x > ^ have particularly 
poor fitting factor This is consistent with the results 
we have seen for the case of non-precessing binaries, as 
shown in Fig. 6, suggesting that the poor performance of 
the non-spinning templates is due to neglecting the sec- 
ular spin-dependent effects in the phase evolution. For 
the case of binary neutron stars, as much as 25% of the 
SNR can be lost if ;^ ^ 0.1. As the mass ratio increases 
(77 «; 0.25) other precessional effects also start to affect 
the performance of the template family. 



VII. CONCLUSIONS 

In this paper, we presented a PN detection template family 
of gravitational waveforms from inspiralling compact binaries 
with non-precessing spins. The waveforms are reparametrized 
in such a way that all spin-dependent terms are described by 
a single "reduced-spin" parameter in an approximate fashion. 
We have shown that the template family has very high overlaps 
with non-precessing binary signals with arbitrary spins and 
mass ratios. The family is also effectual for the detection of 
a significant fraction of precessing binaries in the comparable- 
mass regime. This is due to the fact that, in the comparable- 
mass regime, the secular spin-dependent effects dominate the 
phase evolution (not the oscillatory effects), which are cap- 
tured by the non-precessing templates. Since non-spinning 
templates neglect this effect, a significant fraction of SNR can 
be lost if non-spinning templates are used in the search for 
spinning binaries. 

The simple, closed-form expression of this frequency do- 
main template family makes it very easy to implement this in 
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FIG. 10: Scatter plots of the fitting factor (indicated by the color of the dots) of the reduced spin template family with precessing PN binaries 
as a function of different parameters describing the initial configuration of the binary. The angles and t^r describe the initial orientation of the 
binary with respect to the detector, 11^211 is the spin magnitude of the more massive compact object and 6s2 describes its orientation with respect 
to the orbital angular momentum (see Sec. Ill A for a complete description) while x = {^ - T6r]/ 1 1 3) x.t 'La- + 6Xc, ■ Lai is the initial value of 
the reduced spin parameter. The total mass m/Afg and mass ratio q of the target population in each column are shown in the title. 



search pipelines, and will provide considerable improvement 
over the non-spinning templates. Since a parameter-space 
metric [24] can be easily worked out for this waveform fam- 
ily, the current techniques of template placement and multi- 
detector coincidence tests [21] can be readily used in such a 
search. The parameterization of the inspiral waveforms using 
the reduced spin parameter ;t' can also be extended to the con- 
struction of inspiral-merger-ringdown waveform templates for 
binary black holes, such as the ones presented in [18, 86]. 
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